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1 Introduction 



1.1 Microscopic model 

Each site of a one-dimensional lattice (Z) is initially occupied by one particle A 
with probability p or empty with probability 1 — p. The A particles perform simple 
(continuous time) random walk with a diffusion constant D. We further impose 
the exclusion constraint, namely, each site is occupied by at most one particle. In 
addition the particles may undergo several reaction processes : 

coagulation : A + A —>■ A at a rate k; 
branching : A ^ A + A at a rate A; (1.1) 
decay : ^ — > at a rate 7. 

Owing to the exclusion constraint particles react when they sit on neighboring sites. 
Similarly diffusion takes place only when empty sites allow it to. Coagulation, 
branching and decay, along with diffusive motion, define a reaction-diffusion pro- 
cess that has already received considerable attention in the past : this is the Schldgl 
autocatalytic reaction, which is known to belong to the universality class of directed 
percolation. The d-dimensional generalization of this model has been studied via 
renormalization group techniques by means of e-expansion in the vicinity of the 
upper critical dimension = 4 (|[l]]). The very few analytic results that exist in 
one dimension are based on short-time series expansions ([0]). Our aim in this 
work is to provide a systematic approximation scheme specific to d = 1. 



1.2 Mean field for directed percolation 

It is possible to write a mean field equation for the average local particle density 
n{t) at time t : 

^ = (A-7)n-(A: + A)n2 (1.2) 
From this equation one predicts that in the steady state 

n{oo) = j '[^ > ^ (1.3) 

[ otherwise 

Hence mean field predicts a continuous transition between an active state for 7 < A 
in which a finite fraction of A's survives indefinitely, and an absorbing state in 
which A's have completely disappeared forever, which occurs for 7 > A. In the 
following we shall use 7 as the control parameter and fix all other parameters. At 
the mean-field level we see that the steady state of this system undergoes a second- 
order phase transition between an active state with nonzero A density, at 7 < A, 
and an absorbing ^-free state at 7 > A. Within the mean-field picture the transition 
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occurs at the critical value 7c = A. It is possible to summarize the scaling properties 
of the particle density in a single formula 

n(t) = 6-^^(6-% 6^/^7-7^1) (1.4) 

which holds for 6 » 1 with the arguments of T fixed. This scaling relation defines 
the critical exponents r/, z and v. Their mean-field values are = 0, z = 2 and 
V = 1/2. In the steady state the density behaves as I7 — 7c|^ as 7 — > 7", which 
defines the exponent /? = v(\ + r\)j2. 



1.3 Motivations and outline 

From the analytic point of view there are very few exact or even approximate 
results on directed percolation in low space dimension (see ^ for recent re- 
views). Owing to its ubiquitous nature in the study of stochastic processes, di- 
rected percolation has become the paradigm of out-of-equilibrium systems pos- 
sessing a second-order phase transition in their steady state. Our aim is to remedy 
the scarcity of analytic techniques specific to low and physically relevant space di- 
mensions. Indeed d = 1 is the relevant dimension for instance in the study of sur- 
face growth phenomena. Another motivation comes from particle physics. There 
the branching-coagulation language is used as a phenomenological description of 
hadronic high-energy scattering processes with d = 2 being the physical dimen- 
sion corresponding to the number of transverse space directions. Other applications 
include the study of intermittency (see e.g. Henkel and Peschanski [^), either in 
the Schwinger mechanism or in turbulence [^. Fluctuations play an increas- 
ing role as the dimension is decreased below the upper critical dimension dc = 4, 
which provides further motivation to focus on low dimensions. In the absence of 
any exact solution we believe that our method provides new insight into the pecu- 
liarities of one-dimensional directed percolation. 



The outline of this paper is as follows. In section 2 we map the master equa- 
tion that describes the reaction-diffusion process Eq. (LI) first onto a spin-chain 
problem. The spin-chain is then mapped onto a fermionic field theory. The proce- 
dure, which we describe in great detail, consists in building a fermionic field theory 
starting from a non-hermitian hamiltonian originating from the stochastic process 
Eq. (|ni). This raises a number of difficulties which, to our knowledge, appear for 
the first time in the literature. These are exact mappings. We then proceed with 
the analysis of the full theory describing directed percolation using renormaliza- 
tion group techniques. Our calculations yield numerical estimates for the critical 
exponents. 



2 From the master equation to a quantum spin chain 
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2.1 Master equation 

A microstate of the system described in paragraph is characterized by the 
set of occupation numbers {nj}j^i defined by 

1 if site i is occupied by an A 
if site j is empty 

We now define the spin variable Sj = 2nj — 1. Let s = {sj} denote a generic 
microstate of the system and let it index a set of vectors \s) in a Hilbert space. 
The master equation for the probability of occurence P{s, t) of state s at time t is 
equivalent to an evolution equation for the linear combination 

|$(i))^5]P(s,t)|s) (2.2) 

s 

which reads 

= -mm) (23) 

where H is an evolution operator (also abusively called a Hamiltonian) acting in 
the Hilbert space spun by the states | s) . We introduce the Pauli matrices ai defined 

and also define the raising and lowering operators 

af^]^{af±tay) (2.5) 

We thus have the identities af = 2afal — 1 and ofa^ + (t^ct^ = 1. The variable 
Sj is the eigenvalue of af. One may verify that H may be written in the form 

H — -^diffusion ~l~ -f^decay ~l~ -^^branching ~l~ -^^coagulation (2.6) 

where we have set 

-^diffusion = — — ^i-^i+l (2-7) 

i 



H^...y = lY.^\at-o-) (2.8) 



A 



^branching = ^^(cTf af+i + 2^7+ + a+^jf+i + a+ ^^f ) (2.9) 

i 

k \ 

^coagulation = ^ /Z^'^^i ~ + " " ^i+l^i) (2-10) 

i 

Note that we have dropped all constant terms in Eqs. (2.7- 2.10| ). They ensure the 
conservation of probability but will however play no role in the subsequent analy- 
sis. 
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2.2 General properties of the spin chain and average of observables 

In this paragraph we recall for completeness some of the properties of the spin 



chain Hamiltonian defined by Eq. (P^l). We need to introduce a projection state (p| 
defined by 

(pI^E^^I ^2.11) 

s 

Given a physical observable A{s) we denote by A the operator obtained by replac- 
ing in the explicit expression of A the variable Si by the operator af. For instance 
the choice A{s) = ^{sj + 1), which is the local number of particles at site j, leads 
toi = i(cj| + l).The average of the observable A{s) may be expressed as 

{Aism = {p\A\m) (2.12) 
as was first noticed by Felderhof 

Conservation of probability imposes that (p| is a left eigenstate of H with 
eigenvalue : 

{p\H = (2.13) 

from which it follows that 

yt, (p|i>(t)) = i (2.14) 

Besides, H has at least one right eigenvector with eigenvalue 0, which describes 
the stationary state of the system. The eigenvalues of H all have a positive real 
part. Other details may be found in the reviews by Alcaraz et al. |Q] or Henkel et 



al. QlOp. For our purposes we need one more property of the projection state. It is 



based on the following identity : 
After noticing that 

(p| = (-l|e^i''7 (2.16) 
it becomes possible to express the average of an observable A{s) in the form 

{A{s)){t) = (-l|i'e-^*e^^"7|$(0)) (2.17) 

where A' and H are deduced from A and H as follows. Express these operators 
only in terms of the cr+'s, cr~'s and a~^a~'&. For each j replace in the resulting 
expression by ~'^J ~ '^'^'j^J + 1- This yields the operators A and H. In the 
expression of A one first puts all the 's to the left of the a~ 's then one formally 
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Eq. (2.6) one finds after straightforward manipulations : 



2hi 



sets the cr+'s to 0. This yields the operator A'. With the particular expression of H 

forward 

« = -(" + |)E[' 

i 

+ (7-A) J^ni 

i 

+ {X + k - 2D)'^hihi+i ^2 18) 

i 

i 



X 
2 



where we have adopted the notation hi = <jfa~. As we have already mentioned 
we have omitted a constant term in H. In order to find this constant, one would 
have to push the cr+'s to the left of the d^'s then formally set to all the a'^'s, 
which produces an operator H' and adjust the constant so that {—l\H' = 0. 



3 Fermionic field theory 

3.1 Jordan- Wigner transformation 

We define a set of fermionic creation and annihilation operators cj and Cj by 



where 



4 exp ^i-K 



rii 



■3 ' 



exp 



^7r 



Ci 



(3.1) 



c^-Cj is the particle number operator. It has the eigenvalue 



Uj. This is the Jordan-Wigner transformation. The commutation relations between 



the Ci and cj are 



{ci,c]] = 6ij, {4,c]} = = {ci,Cj} 
An equivalent expression for the string operator in Eq. ( [3^] ) is 

As an example, for diffusion alone, the evolution operator would read 

-^diffusion = -D ^ [cjcj+i + cf^iQ - 2hi + 2hihi+i 



(3.2) 
(3.3) 

(3.4) 



It is surprisingly left invariant by the transformation Eq. (2.15). In the quadratic 
part of Eq. (3^) one may recognize the diffusion process, that results in a Lapla- 
cian. However diffusion takes place with the constraint of excluded volume, which 
accounts for the quartic interaction term. 
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3.2 Coherent states 

As we have seen in the previous section, the average of a physical observable 
can be cast in the form 

{A){t) = (-l|i'e-^*eS>'^>"|$(0)) (3.5) 

We shall restrict our analysis to the case of an initial state with particles randomly 
and independently distributed with local density pj at site j, that is with the initial 
distribution 



(3.6) 



P(n, t = 0)=Y[ [pjSn^,l + (1 - Pj)Sn,,0 

j 

Under those assumptions it is easy to check that 

|$'(0)) = eS>'"»^|^>(0)) = e^^-^^'4| - 1) (3.7) 

At this stage we will denote by |0) the vacuum (i.e. particle-free) state, instead of 
I — 1), as the former notation is more appropriate in the particle number language. 
To summarize, the average of an observable A reads 

{A){t) = (0|i'e-f^*|$'(0)) (3.8) 
We rewrite the exponential factor in Eq. (|3.8[) using Trotter formula : 



e-^*= lim (l-$^Y (3.9) 

We now introduce the coherent states associated to the creation and annihilation 
operators. Let \ipij), for j = 0, and i G Z denote a set of coherent states 
associated to q and c| : 

{iPijl = {0\e-h'^'j'^'^e-^''^'j (3.10) 

These coherent states are indexed by a pair of conjugate Grassmann variables 
and By definition they have the property that 

Ci\ipij) = ipijli^ij) , {Aj\4 = i^ijl'^ij (3-11) 



With the definition Eq. ( |3.10| ) one may verify that 

{iJij\^P^J) = 1 (3.12) 
and that, with the integration convention f dip ip = f dip* ip* = 1, one has 



/ 



dV'*.d^,,e-'^?.^-C,V'r = ^nn. (3.13) 
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which impUes 



Eq. ( [3.14] ) generaUzed to the state 



(3.14) 



(3.15) 



leads us to 



-Ht 



N 



j=0 i 



N-\i 



(3.16) 



(V'i|i-^^|V'o)(V'ol 



Thus we need the quantities 



as well as 



(3.17) 



(3.18) 



In order to obtain Hjj^i one replaces in the normal ordered form of H the annihi- 
lators C£ by i^ej-i and the creators Cm by ip^j. The projection onto the initial state 
|$'(0)) reads ' 

('0i=o|^'(O)) = ...e''^'^«e'''+i'^?+i.o... = Pi['il;^{x,0)] (3.19) 
In a similar way, one finds 

(0|i'|V',=jv) = A'{%bN) (3.20) 
where A' is some function of V'at the precise determination of which we postpone 



to paragraph 3.5 



3.3 Continuous field theory 

It is convenient to split Hj^i j in the form 



(3.21) 



where ^j+ij and Oj+ij denote the even and odd components of Hj^ij, respec- 
tively. In the limit At — > the ordered product 



tr = (1 - AtHN,N-i).:il - At/7i,, 



(3.22) 
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does not reduce to a simple exponential as the terms in the product do not commute. 
We may use an identity ([[TTI]) for a sequence of odd Grassmann variables ai,..., 

OtN ■ 



,E„ an+E^=l Em = l """" 



Hence we may rewrite tl in the form 

tr = exp 



E ^^+1.^ + E E 

j j=o e=o 



'e+i,i 



(3.23) 



(3.24) 



which yields in the continuum limit : 



: exp 



+ 



dtH[m,rit)] 



dMi'e(t-t')O[^(t),V*(t)]O['0(i')>^*(O] 



(3.25) 



We define an action S by 



dx dt ijjdtilj 
dtH[^p,4)] 

dt dt' Q{t - t')0[ij{t),i){t)]0[i)(t'),i){t')] 



(3.26) 



where the notation '^(t)] means the (normal-ordered) H in which the cre- 

ators and annihilators have been replaced by their Grassmann eigenvalues after the 
former have been moved to the left of the latter. The quantity 0[V'(t), '^(i)] is the 
odd component of H[^{t),'il){t)]. Thus we come up with an action that comprises 
terms that are nonlocal in space and time. In what follows it could be possible to 
proceed entirely within the fermion operator formalism (pertubation expansions 
have been checked this way). There one would have to use time-ordered products 
in evaluating correlation functions, the effect of which is precisely to rebuild effec- 
tively nonlocal interactions. 



Again using Eq. ( [3.23| ), we find 

Pi[V^(x,0)] = exp [Y^pMQ) +Y,PiP3^ii^)^M 



(3.27) 



l<3 



which allows us to express the average of the physical observable ^ in a path 
integral form as follows : 

{A){tf) = j V4,VijA'Vl'{x,tf)]&-^^^^^'^Pi[4,{xM (3.28) 

We are now in a position to apply the existing techniques of field theory. 
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3.4 Comparison with bosonic field theories for stochastic processes 



To conclude this paragraph we would like to make a comparison with reaction- 
diffusion problems in which the exclusion constraint is not imposed ([|T2[ |l! 



There bosonic operators (aj, a|), rather than spin operators, are used to express the 
evolution operator H. The projection state reads 

(p| = (0|e^»''' (3.29) 



and the bosonic counterpart to identity Eq. ( 2.15 ) is e'^'a^ = {a\ 
commutation of the factor e^*"* through the e^^* can certainly be performed 
from the outset, but is not compulsory. It is possible by an appropriate shift of the 
field associated to the creation operators to perform this transformation directly 
within the functional integral formulation. In the present case it seems that such 
an a posteriori nonlinear change of Grassmann field having the same effect as the 



nontrivial transformation Eq. (2.15) would be difficult to find and to handle. 



3.5 Expectation values of some observables 

The main observable of interest is the local particle number A = rii. Hence 
A = hi and A = hi + a~, so that A' = . Since 

(0|i' = {0\ci (3.30) 

we can write the average of the local particle number at site i as the functional 
integral 

{ni)it) = j PV^PV ^i{t) e-^[^'^lPi[V^] (3.31) 

Another physical quantity of interest is the equal time two-point correlation 
function A = niUj, with i < j. The associated A' reads 

A' = a-aj = JJ(1 - 2n^)cjQ (3.32) 

hence (OjA' = (0|cjQ which yields the equal time correlation function 

{ninj){t) = j Vi>V^ ^j{t)tlJi{t) e-^t'^'^^lPi (3.33) 

For i = j the above formula does not hold, and one finds A' = Cj so that one re- 
covers (n|) = {rii). Other equal time averages can be derived in a similar fashion. 
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4 Directed Percolation 



4.1 Action 

We now return to our problem of interest, directed percolation. It is described 



by the evolution operator Eq. ( |2.18[ ), which splits into several contributions. We 



find it convenient to introduce the non local functional 

^^(f) = J|e-2V;K<)^£W (4.1) 



which takes into account the string operator present in the definition of the Jordan- 
Wigner fermions. The action describing the whole reaction-diffusion process reads 



Sdp = y dt i;i(dt + ^-X-{D + ^)A^ 



+ (X + k- 2D)il)i'ipiil)iJ^iil)i+i 

- (A + k.)4)i^iii-idtpi 

- ^ / mt'e{t - 1') [-(A + k)i/jii/ji^i-idi^i + xi^ii^i^i-idiJi] (t) x 

[-(A + k)i)jtpjij-idiijj + \i)j^jij-idi)j\ {t') 

(4.2) 

with dfi = ^{fi+i - /j-i) and A/j = /j+i - 2/j + Finally note that we 
have omitted terms (located at the time slice t = 0) describing the initial state of 
the system. 

4.2 Scaling analysis 

After rescaling fields and coupling constants we may take 



S[^, iP]= j dxdt [ijj{dt + Da- Ddl))^ + DgiP^PCid^iP - d^ijj)] 

+ / dxdydtdt' [ui{i:^jjCdx'4'}{x,t)Q{t - t'){i>'>pCdyi)}(y,t') 

_ _ , _ (4.3) 

+ U3{'ilj'ip£,dxij}(x,t)&it - t'){i^'iPidyip}{y,t') 

+Ui{'4'i}idx^]i^^^t)Q{t - t'){'4}i)idyi)]i^y^t')\ 

as the starting point of the subsequent analysis. The coupling constants are ui = 
U2 = —U3 = —U4 = {DgY- We find the field and coupling constants dimensions 
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by power counting in momentum unit : 

m = m = \. [x^] = [t] = -2, [g^] = l (4.4) 

There is an interesting difference with the one-dimensional Reggeon field theory : 
there the equivalent coupling constant has dimension [g"^] = e = 4 — d = 3. It 
is natural to expect, then, that perturbation will lead to more sensible results than 
the mere d ^ 1 extrapolation of the e-expansion. 



In the action Eq. ( |4.3[ ) the nonlocal exponential factors 

=exp(-2 / dyV^V(y,i)) (4-5) 

J — oo 

play a crucial role. Indeed, by construction their presence breaks the space x ^ —x 
symmetry of the action, along with translation invariance. This signals that the right 
quantities to study in a renormalization perpective are not the usual vertex func- 
tions. 



Finally note the absence of the quartic interaction term of the form 

j dt'Y^ipiipi'ilji+i'ipi+i (4.6) 

i 

It may readily be seen that this term can be written in the form 

ydt J^V^.V'.V^mV'm =/di|r^S^ 

i 

X {2TT)6{qi + 92 + 93 + (?4)(e'(''^+''^) - 1 - i{q3 + qi)) 
X V'(9i,i)V'(92,i)^fe,i)V'(94,t) 

(4.7) 

which, in the small momentum region yields an irrelevant operator of negative 
dimension equal to -1. 



4.3 Correlation functions and renormalization 

At this stage we have to choose a renormalization scheme. We have explained 
why the natural choice of examining the divergences of the vertex functions is not 
adapted here since we are dealing with a theory that neither has translation invari- 
ance nor space reflection symmetry. This is due to the anisotropic construction of 
the fermion operators. In fact the two-point vertex function can still be used to 
renormalize the fields, the diffusion constant and the mass. As for the coupling 
constant we have chosen a particular class of correlation functions in which space 
arguments are carefully ordered. As could be expected the perturbation series or- 
ganizes in powers of the dimensionless coupling g'^ / y/a, and our interest lies in 
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the critical regime cj — > 0. Renormalization is thus needed to extract physical in- 
formation from the above naive expansion. 

We will renormalize the massive theory at zero external frequency and momen- 
tum ([14 1). We define the renormalized fields and parameters by 

D = Z-^ZdD^ 

g = Z^^Z-^/^Zgg^J^ (4.8) 
The first three Z-factors are determined by the following conditions : 

^^f-y(0,0) = l, -^^(0,0) =D^, rl^'^> (0,0) =D^a^ (4.9) 

The one-loop expression for the propagator, calculated after the diagram depicted 
in Fig. [ij reads 



+ Da + g ' 



2TTq^ + a 

An explicit computation of the integrals yields 



+ Da + : 



2 f dg 

27r g2 + o- 



(4.11) 



In order to renormalize the composite operator ipip one notes, as in [15] that a 
vertex function with a ipil) insertion is obtained by differentiating the corresponding 
vertex function with respect to Da and we choose the condition r[^'^'^^ (0, 0) = 1. 
As for the renormalization of the coupling constant g we may choose to look at the 
correlation function 

G(^'2)(rci,ti;x2,t2,X3,t3) = {^{xiM)i^{x2M)'^{xi,h)) (4.12) 

with the restriction that X2 < xi < x^. Following the physical explanation pro- 
posed by Cardy [15], one understands it is more convenient to look at the the space 
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integral of G^^'^^ with the space arguments kept ordered as specified. We now 
define a quantity very similar to a vertex function : we first continue the latter 
function correlation function in which the arguments are ordered to unconstrained 
space arguments. This does allow for an easy calculation of its Fourier transform. 
We define the associated IPI function as r(^'^)(A;i, A;2, ^3) (by convention, leg 3 
carries the derivative ik^). We emphasize that this is not the usual vertex function 
as defined by the Legendre transform of the logarithm of the partition function. We 
define the renormalized coupling g^^ by 



({0,0}) =2L»r5r 



(4.13) 



d{ih) 

We refer the reader to Fig. ^ for the corresponding one-loop diagram. A rather 
tedious calculation yields 

vr (g2 + aY 4 



diih) 



({0,0}) =21)5 1 



2 



dq q^ 

2^(g2+a)2 



2Dg 



1 



(4.14) 



so that the Z-factors read 



1 + 4, Zd 



1 + — , 
16 



Za 



Denoting the renormalized correlation length by ^ 
functions 7^ = i = ^, D,a, g follow : 



1 + ^ (4.15) 



— 1/2 

ctr , the one-loop Wilson 



The /3-function reads 
hence, to one-loop. 



9l 



9l 



9l 



16 ' ~ 4 ' '^9 



9k 



t (4-16) 



(4.17) 



-9k 



3 2 
5r 



32 



(4.18) 



The renormalization group flow then has a single stable fixed point g'^ = The 
critical exponents can be expressed in terms of the Wilson functions at the fixed 
point : 

2 



2-7B + 7; 



3 

'^"^ = 2 - 7d + 7: ^ ^ 



(4.19) 
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The above numerical values are in reasonable agreement with the precise numerical 
estimates of Dickman |jl9l] (77 = -0.405, z = 1.767, 1/ = 1.142). 



5 Conclusions 

5.1 Comparison with existing approaches 

Of course the method we have presented suffers two obvious drawbacks. First 
the numerical values of the critical exponents poorly compare to those given by 
numerical simulations. This is due to the absence of a small parameter validat- 
ing the first order loop expansion. From that point of view the situation is quite 
similar to the low-dimensional extrapolation of e-expansions. A related question 
concerns the convergence of the series expansion of the exponents in terms of the 
fixed-point renormalized coupling. Second we lack a decisive argument in favor of 
the renormalizability of the theory, a property we trust in, but that we were unable 
to establish on a rigorous basis. Arguments in favor of the existence of consistent 
renormalization rely on the locality of the original spin theory. The scheme that 
we have introduced, based upon the study of space ordered correlation functions 
(and their related loop integrals), must now be extended to all orders. However we 
emphasize that we have built an entirely new and systematic analytic technique 
that provides to our knowledge the only alternative to e-expansion of the Reggeon 
field theory and to short time series expansions. The approximation scheme we 
have developped has in principle a wide range of applicability and exhibits many 
analogies with its bosonic counterpart. Its main advantage is that it naturally in- 
corporates the simplification induced by the one-dimensional topology. As a final 
comment, we find it amazing enough that the present field theory has appeared in 



a slightly different form in a quite different context [15]. 



5.2 Prospects 



We have also applied the above formalism ([16]) to the much studied annihi 



lation reaction ^ + ^ — > 0. The methods of integrable spin chains ([y, |10|]) work 



at a particular diffusion contant-annihilation rate ratio (which corresponds to a free 



fermion Hamiltonian Eq. (2.6)) and leads, as expected, to a l/yi decay of the par- 
ticle density. The present method very easily allows to extend this result to any 
value of the parameters. 

There exists another class of reaction-diffusion problems, namely branching 
and annihilating random walks ([pT]]), in which the formalism we have developped 
should profitably apply. We expect it to provide a better qualitative picture in one 
dimension than previous attempts ([|l7]]) using bosonic theory down from the up- 
per critical dimension. It also remains to see how the formalism can be extended to 



cope with several species reaction-diffusion problems (e.g. []18|]) for which no an 



alytical nor numerical results exist in one dimension. Our tool should prove useful 
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in the study of those problems, where other techniques have, up to now, failed. 
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Fig. 1: A plain leg denotes a tp field while an arrowed one denotes a conjugate ^. 
This is a diagram contributing to the renormalization of the propagator. Two of the 
four internal legs carry a derivative. In Fourier space, the loop integral has to take 
into account the product of the momenta carried by those legs. 



FIGURES 
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Fig. 2: One of the diagrams entering the renormaUzation of g. It involves ui and 
9- 
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